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Abstract This review describes the observational properties of radio pul- 
sars, fast rotating neutron stars, emitting radio waves. After the introduc- 
tion we give a list of milestones in pulsar research. The following chapters 
concentrate on pulsar morphology: the characteristic pulsar parameters such 
as pulse shape, pulsar spectrum, polarization and time dependence. We give 
information on the evolution of pulsars with frequency since this has a di- 
rect connection with the emission heights, as postulated in the radius to 
frequency mapping (RFM) concept. We deal successively with the proper- 
ties of normal (slow) pulsars and of millisecond (fast - recycled) pulsars. 
The final chapters give the distribution characteristics of the presently cat- 
alogued 1300 objects. 

Key words Radio pulsar morphology: pulse shapes - spectra - polariza- 
tion - distributions 



1 Introduction 

The first publication [54] announcing the discovery of a new class of ob- 
jects, soon to be known as pulsars [99] [90], appeared in the literature al- 
most exactly 35 years ago. Since then, large strides toward the understand- 
ing of pulsars and their radiation have been made. We know that they are 
fast rotating neutron stars [3] [4] [118] with masses of the order of 1 Mq 
[122] [29] [151] [58] [20] [149] [144] [49] [57] and strong magnetic field (between 
10^ and 10" Gauss [39] [15] [138] [109] [2] [169]. They are probably born dur- 
ing a Supernova Type II explosion of a (massive) late-type star (Fig. 1). 
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Fig. 1 Scenario of the birth and evolution of pulsars 



Pulsars emit highly accurate periodic signals (mostly in radio waves), 
beamed in a cone of radiation, centred around their magnetic axis. These 
signals reveal the period of rotation of the neutron star, which radiates, 
like a light-house, once per revolution. The lighthouse effect is caused by 
the dipolar magnetic field not being aligned with the rotation axis of the 
neutron star. As a consequence of the magnetic field, pulsar radiation is 
highly polarized. Their period of rotation (P) varies between 1.557 ms (642 
Hz) [8] and 8.5 s (0.12 Hz) [168]. As pulsars rotate, they loose energy and 
their rate of rotation decreases. This period derivative (P) is an important 
observational parameter. In 1975 a pulsar in a binary system was discovered 
[56]. Slow rotating pulsars, with rotation period P > 20 ms and P > 10~^* 
are considered to be normal pulsars. Pulsars with P < 20 ms and P < 10~^* 
are called millisecond or recycled pulsars. A full list of 1300 pulsars [105] is 
available at http://www.atnf.csiro.au/rcscarch/pulsar/psrcat/. 

The period of rotation of normal pulsars increases with time, an obser- 
vational fact, discovered during the very early history of pulsars [134] and 
led to the rejection of suggestions that the periodic signals could be due 
to the orbital period of binary stars. The orbital period of an isolated bi- 
nary system decreases as it looses energy, whereas the period of a rotating 
body increases as it looses energy. Millisecond pulsars are considered to be 
recycled pulsars, spun up by mass transfer (accretion) from a binary com- 
panion [1]. In the early history of pulsars, models involving pulsations of 
white dwarfs and neutron stars were also proposed and quickly rejected. 
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The discovery of the Vela pulsar (PSR 0833-45 - see Table 1) led to the 
suggestion of pulsar - supernova association. This suggestion was corrobo- 
rated by the discovery of a pulsar (PSR 053H-21) in the heart of the Crab 
supernova remnant with a period of 33 milliseconds (see Table 1), which led 
to the unequivocal association of (radio) pulsars with rotating neutron stars 
[45]. Subsequent polarimetric observations led to the establishment of the 
"Rotating Vector Model" [125]. Soon after its discovery, the Crab pulsar 
was post-detected in earlier (1967, Jun 04) archived X-ray data [37] and 
soft 7-ray data [38]. To date (excluding non-pulsed detections of neutron 
stars) 5 normal pulsars have been detected in the optical, 17 normal and 
6 millisecond pulsars in X-rays and 7 normal pulsars in 7-rays (up to 10^^ 
Hz, covering, thus, the largest frequency range of all known compact species 
emitting intense radiation in the Universe) (updated information from [13] - 
W. Becker, private communication). The first (and fastest, up to now) mil- 



Table 1 Pulsar milestones 



Date 


Milestone 


Reference 


193z, i^eb 


The discovery of neutron 


Fool fool 

[22], [23] 


1934, 1967 


Neutron stars are predicted 


[3], [121] 


1939 


Neutron stars Equation of State 


[119] 


1967, Nov 28 


The discovery of pulsars 


[54] 


1967, Nov 28-1968, Mar 03 


"Pulsar" designation 


[14] 


1968, Mar 15 


First published "Pulsar" designation 


[154], [152] 


1968 


Discovery of the Vela pulsar 


[81] 


1968 


Discovery of the Crab pulsar 


[145], [25] 


1968, Feb 24 


Dispersion measure measured 


[30] 


1968, Apr 01 


First polarimetric observation 


[89] 


1968, Apr 03 


Faraday Rotation measured 


[143] 


1968 


Gravitational emission proposed 


[158] 


1968 


Lighthouse model proposed 


[45] 


1968, Nov 


Galactic distribution established 


[82] 


1969 


Post-detection of pulsar X-rays 


[37] 


1969 


Post-detection of pulsar 7-rays 


[38] 


1969 


Scintillation explained 


[135] 


1969 


Rotating Vector Model proposed 


[125] 


1969 


Observations of pulsar glitches 


[126] 


1969, Aug 


First emission process proposed 


[46] 


1974, Jul 02 


Discovery of binary pulsars 


[56] 


1975 


First complete theory attempted 


[139] 


1982 


Discovery of millisecond pulsars 


[8] 


1991, Sep 15 


Detection of extrasolar planets 


[163] 


1992, Oct 19 


Detection at mm-wavelengths 


[162] 


1998 


Discovery of magnetars 


[69] 


1998, Nov 05 


Discovery of the lOOO' pulsar 


[28] 
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lisecond pulsar, PSR B1937+214 was discovered in 1982 ([8]. Ten years later 
the first planetary system (two planets orbiting PSR B1257+12), outside 
the solar system, was discovered ([163]). 

The discovery of millisecond and binary pulsars (see Table 1) gave new 
insight in pulsar research. It was soon realized that the extremely fast ro- 
tation of millisecond pulsars could only be explained by transfer of angular 
momentum from companion stars. The notion of recycling of old and ex- 
hausted neutron stars became popular. On the other hand binary pulsars 
were used to check the effects of General Relativity, which successfully sur- 
vived the new tests. For these particular discoveries the Nobel Price in 
Physics (see http://www.nohel.se/physics/laureates/) was twice awarded to 
pulsar researchers: In 1974 it was awarded to Antony Hewish for the dis- 
covery of the first pulsar, PRS B1919-h21, and in 1993 to Russell A. Hulse 
and Joseph H. Taylor for the discovery and subsequent work on the physics 
of the binary pulsar PSR B1913-I-16. Recently the discovery of a double 
pulsar system, PSR J0737-3039B (spin period 2.7 sec), PSR J0737-3039A 
(spin period 22 msec) in a 2.4 hr eccentric orbit was announced [86]. This 
system will, certainly play a very important role in deciphering the pulsar 
riddle and in testing physical theories. 

Meanwhile, a large number of strange and unexpected properties of pul- 
sar emission were observed. Drifting subpulses, mode changing and nulling 
were among the first such properties to be studied[129]. Extremely narrow 
pulses [51] [76] [52] were soon to become an important tool for theoretical 
investigations and 90° polarization jumps [orthogonal modes) imposed re- 
strictions to the existing models. In addition, period glitches were observed, 
during which the pulsar period decreased by a large amount and then, within 
a few days, it increased again to its previous value. Glitch properties are 
used to study the physical structure of neutron stars. 



2 Pulsar Milestones 

Before proceeding with the description of the characteristics of pulsar emis- 
sion, it is worth looking back and paying tribute to the main discoveries 
concerning pulsars. Table 1, summarises the most important steps in pulsar 
research, both in observation and theory. 

It is not always easy to unearth a "first" date or a "first" publication. For 
example, J. Bell-Burnell has communicated to us her personal view [154] 
and we have found a report in a March 1968 issue of the Time magazine 
[152] referring to the first use of the word pulsar. Nevertheless, we believe 
that still earlier references to the "pulsar" notation may have escaped our 
search. The list can be expanded to include several other "firsts" . However, 
we decided to restrict its size to what we consider to be the most important 
ones. 
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3 Morphology 

Pulsars present to the observer a most complicated set of time variable 
phenomena. On the one hand, their period of rotation, P, varies from 1.5 
millisecond up to 8.5 seconds from object to object. On the other, this period 
increases regularly with time, as the pulsar looses gravitational energy and 
slows down. The deceleration is expressed quantitatively by the measured 
P. The youngest (normal) pulsars exhibit usually the largest deceleration 
and thus they demonstrate the largest P. Quite to the opposite, millisecond 
pulsars have very low P and are interpreted to be recycled objects, slow, 
normal pulsars that have been sped up by accretion from a binary coun- 
terpart. The deceleration process is not always constant. Some pulsars are 
known to exhibit glitches, sudden acceleration to higher periods, which are 
interpreted to be due to crust related "starquakcs" . 

Taking into account single pulse intensities and their duration, extremely 
high brightness temperatures (of the order of 10^^ K) are calculated, espe- 
cially for their low frequency emission. This constrains their radio emission 
to be coherent (see also [83] and references therein). The intensity of sin- 
gle pulses show enormous intrinsic variations. The situation is even further 
complicated by the fact that interstellar scintillation introduces an addi- 
tional fluctuating effect, in particular at lower radio frequencies. However 
with enough repeated observations the average emitted pulsar flux density 
can be determined. The addition of a sequence of single pulses leads quickly 
(usually, after a few hundred pulses) to a stable pulse shape, a signature of 
the geometry of lighthouse emission mechanism. The pulse shape is a firm 
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Fig. 2 A sequence of 17 single pulses observed simultaneously at 4 frequencies, 
plotted one after the other. Courtesy of A. Karastergiou 
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characteristic of a pulsar, having anything from one up to nine sub-pulse 
components. Within the sub-pulse structure a number of phenomena has 
been observed. Sub-pulse drifting (see above) reveals the emitting beam 
structure. In addition, the existence of very narrow pulsar micro- and nano- 
structure [35] [52], that are the signs of individual emission regions, has been 
confirmed. Pulsar radiation is highly polarized in a most complicated way. 
At low radio frequencies some pulsars are almost 100% linearly polarized. 
Others have very high and variable circular polarization. The development 
of polarization with frequency is radically different from all other radio 
sources. The polarization may be high at low frequencies while dropping 
rapidly to zero at high frequencies. Possibly this is a hint for a coherent 
(low frequencies) - incoherent (high frequencies) emission mechanism, an 
effect corroborated by high frequency pulsar spectra [72] . 




0.0 3.0 S.l g.l 12.1 15.2 DEG 

Fig. 3 (Top): A sequence of 100 single pulses from PSR 1133+16 plotted un- 
derneath each other (Noise wings are suppressed) . (Bottom): By adding the above 
single pulses, we get the Integrated profile 
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3.1 Displaying pulses 

Pulsars are immediately recognised from the periodic nature of their radia- 
tion. Pulse sequences (Fig. 2) can be depicted in a much more compact and 
informative way if their period is accurately known. Then, instead of show- 
ing the pulses one-after-the-other, they can be displayed one-underncath- 
the-other. Thus, many more pulses can be conveniently accommodated in 
a single graph. By ignoring the unpulsed noise either side of the pulses, 
high time resolution single pulses are readily displayed (Fig. 3 - top). It is 
evident from Figures 2 and 3 that individual, single pulses vary greatly in 
intensity. Most of these variations are intrinsic. Some are due to interstellar 
scintillation. However, if a large number of single pulses are added together, 
a very stable profile is obtained (Fig. 3 - bottom). For most pulsars, this 
integrated profile characterizes uniquely a pulsar at a particular frequency. 
The stability of pulsar integrated profiles has been thoroughly investigated 
for a large sample of normal pulsars [53] and some millisecond pulsars [65] . 



3.2 Integrated Pulses 

In the integrated pulsar profiles distinct components can be identified [140] [67] . 
They are thought to represent coherent physical regions in the magneto- 
sphere of the star. Therefore their properties arc of extreme importance 




Fig. 4 Time-aligned profiles at three frequencies, with gauss-fitted components. 
Radius to Frequency Mapping (RFM) is obvious in this Figure. Left panel: PSR 
B0355-t-54. Right panel: PSR B0740-28. 
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for understanding the emission mechanism of pulsars. These components 
are often blended and their shape and longitudinal location within the 
pulsar profile is difficult to establish. A large collection of pulse profiles 
can be found in the European Pulsar Network data archive at www.mpifr- 
bonn.mpd.de/div/pulsar/data/. Although there are no rigorous theoretical 
arguments, usually they can be fitted with gaussian curves, the parame- 
ters of which are easily obtained. Experience has shown that pulsar profiles 
can indeed be reconstructed by a sum of individual gaussian components. 
This method usually involves some assumptions which can be minimized by 
reducing the number of degrees of freedom of the gauss-fitting procedure 
[71] [70]. Some gauss-fitted components are shown in Fig. 4. 




Fig. 5 Integrated profiles of pulsars at I.4I GHz [I4O] 



Morphology and Characteristics of Radio Pulsars 



9 




-25 -15 -5 5 15 25 (Deg) -30 -10 10 30 (Deg) -50 -30 -10 10 30 (Deg) 




-40 -20 20 (Deg) -20 -10 10 (Deg) -30 -20 -10 10 20 (Deg) 




-25 -15 -5 5 15 (Deg) -20 -ID 10 (Deg) -30 -20 -10 10 20 (Deg) 



Fig. 6 Integrated profiles of pulsars at 10.7 GHz [14-0] 



3.2.1 Normal Pulsars (P > 20 ms) Pulse profiles come in a variety of 
shapes (Fig. 5, 6). In most cases they can be represented by a smooth 
curve with a single (almost gaussian) component or with two or more com- 
ponents. Soon after the discovery of pulsars, it was realised that their in- 
tegrated profiles exhibit important morphological differences. In order to 
explain double profiles, the hollow cone model was proposed in the early 
seventies [68] [6] [120]. Triple profiles were explained by the introduction of 
a central pencil beam and five-component profiles were interpreted by as- 
suming a more complicated beam, comprising of a central beam surrounded 
by an inner and an outer cone [127] [128] [130] [131]. The number of cones 
that can be accommodated within the narrow polar cap region, whose ra- 
dius is bound by the last open magnetic field lines cannot be very large. 
It is interesting that pulsars exhibiting a single core component, seem to 
occupy a distinct region in the PP diagram [43] (Fig. 21). On the other 
hand a patchy beam model was proposed [87], according to which pulsar 
beams are patchy, with components randomly located within the last open 
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Fig. 7 Integrated profiles of millisecond 
tions, M. Kramer et al. 



ars at I.4I GHz. Effelsberg observa- 



magnetic field lines. This model is based on the fact that single pulses vary 
in intensity and often they seem to be missing altogether. There have been 
several attempts to explain the pulsar beam shapes using both theoretical 
and geometrical arguments [44][98][42][110]. However, there are still many 
uncertainties due to the erratic behaviour of pulsar emission and the lack 
of an accepted model for pulsar radiation. 



3.2.2 Millisecond (Recycled) Pulsars (P < 20 ms) The first millisecond 
pulsars to be discovered exhibited rather simple profiles. Nowadays about 
100 millisecond pulsars have been detected, many of which have complex 
integrated profiles [75], not dissimilar to the profiles of normal pulsars (Fig. 
7). One difference is that millisecond pulsars tend to have wider profiles 
than normal pulsars. 



3.2.3 Radius to Frequency Mapping The frequency development of pulse 
shapes has led to the concept of a Radius to Frequency Mapping (RFM 
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Fig. 8 Radius to Frequency Mapping (RFM). High frequeneies are emitted close 
to the surface of the star. Low frequencies are emitted higher up, where the cone 
of emission is wider 

- Fig. 8), according to which, higher radio frequencies are emitted in the 
lower reaches of the pulsar magnetosphere (closer to the surface of the star) . 
In order to investigate the RFM concept, time-aligned pulse shapes are 
needed, that require accurate pulsar timing. Following the standard pulsar 
model [139] the integrated pulse width is expected to decrease monotonically 
with frequency (RFM effect). This effect was implicit in earlier work [142] 
and has been extensively investigated ever since [26] [157]. Multi frequency 
observations of pulsars have confirmed the narrowing of pulse profiles with 
frequency [79] [75]. The effect is adequately demonstrated for both normal 
and millisecond pulsars in Figures 4, 9 and 10. 

3.2.4 The Crab pulsar from radio frequencies to j-rays The Crab pulsar 
(PSR B0531+21) is probably, the best studied pulsar. Soon after its discov- 
ery, archival searches led to post detection of pulsed emission at X-rays and 
7-rays [37] , [38] . Its "main-pulse - interpulse" integrated profile is unmistak- 
ably evident throughout the electromagnetic spectrum (Fig. 11). However, 
carefully time-aligned profiles [111] reveal slight, but significant, displace- 
ment of its high frequency components from its lower (radio) components. 
This has been interpreted as evidence of two different mechanisms of emis- 
sion, with the high frequency emission (optical. X-rays, 7-rays) originating 
in a region close to the light cylinder. Furthermore, recent investigations 
[111] [63] have revealed that between 4.7 GHz and 8.4 GHz extra compo- 
nents appear in its integrated profile (Fig. 11). These components impose 
additional difficulties in the investigation of the emission mechanism of this 
interesting object. 
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Fig. 9 The development of pulsar integrated profiles with frequency. Note that new 
components can appear both at high and low frequencies. Courtesy M. Kramer 



3.3 Single Pulses 

Investigations of single pulses is of utmost importance, showing the vari- 
ability and spatial structure of the pulsar emission process. Already in the 
discovery paper [54] and other early papers [30] [82] [27] [32] [137] single pulses 
at low frequencies were studied. Single pulses show a variety of sharp emis- 
sion structures from millisecond through microsecond down to nanosecond 
range [27] [51] [80] [52]. Early two- frequency simultaneous observations sug- 
gested that the emission is inherently broad-band [9], i.e. emission is corre- 
lated over a wide frequency range. Later observations [16] [12] [31] [64] were 
made for total intensity only and at most for two simultaneous frequencies. 
Most of these investigations were then used to determine the bandwidth of 
the emission. More recently, using many radio telescopes at different fre- 
quencies, simultaneously, at up to five frequencies [77] have investigated the 
single pulse characteristics in detail (see also Fig. 2). 



3.4 Millistructure, Micro structure, Nanostructure 



The fact that very short time structures are present in single pulses was 
noted from the very beginning [54] leading to a stringent requirement for 
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Fig. 10 The development of millisecond pulsars integrated profiles with frequency. 
Adapted from [75] 



theories attempting to explain the origin of pulsar emission. First observa- 
tion of pulsar millistructure was made soon after their discovery [27] . The 
limitation was due to signal to noise in the narrow band receiver needed 
to show such short time structures. A few years later a de-dispersion tech- 
nique at the frequency of 115.5 MHz revealed time structures as short as 8 
US in PSR 0950-h08 [51]. Direct observations at 1420 MHz on PSR 1133-hl6 
[36] resolved structures with time scale of 14 fj,s. The microstructure was 
found to be broad-band [136] [16]. Periodic structures were observed that 
seem to be also correlated across a wide frequency range [17]. More recently 
[80] microstructure investigations were extended up to 4.85 GHz, showing 
that many pulsar have this emission signature. This result was confirmed 
in the studies of the Vela pulsar which shows microstructure [76] in most 
of the pulses. Most recently the time structure studies were taken in to the 
nanosecond range with observations of giant pulses from the Crab pulsar 
[52]. Their best time resolution was in fact 2 nanoseconds. This latest ob- 
servation suggests that the plasma responsible for such emission must be of 
the order of one meter in size. If the emission is isotropic, these nanosecond 
pulses must be the brightest transient source in the radio sky. 
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Fig. 11 The integrated profile of the Crab pulsar from low frequency radio waves 
to soft gamma-rays. Note the absence of the main pulse and the peculiar extra 
components between J^.l and 8.4 GHz. Courtesy T.H. Hankins 

4 Flux densities 

The flux density of a radio source and its frequency evolution (its spectrum) 
are basic information that relate to the emission mechanism. However, one 
of the problems with flux densities is that pulsars vary on various time 
scales: (a) due to inherent variations [147] or (b) due to scintillations [94] 
or (c) due to scattering [84]. 

For low frequency radio waves of the Milky Way the spectrum of the 
radio emission could be explained only by the non-thermal (synchrotron) 
emission process. From the very beginning of pulsar observations it was clear 
that the spectra of pulsars were very different from all other known radio 
sources [88] [137]. Instead of values of a - 0.8 ( S = v""), as in cosmic radio 
sources, the observed spectral index of pulsars was a ~ - 1.5 on average. The 
spectral index at the highest frequencies was for some objects much higher 
than this average value. At lower radio frequencies a spectral turn-over was 
observed. It was known in the earliest papers that pulsars were very weak 
at higher frequencies posing in fact a great instrumental challenge to study 
these objects at cm wavelengths. 
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12 Spectra of pulsars detected at 9 mm. [160] 



4-.1 Normal Pulsars (P > 20 ms) 



The first spectra of pulsars, using flux density values at three (plus an upper 
limit at 1.4 GHz) frequencies, were obtained in 1968 [88]. Earlier observa- 
tions [137], published slightly later, were obtained at five frequencies, four 
of them simultaneously, giving an average spectrum and spectra of individ- 
ual (single) pulses of PSR1919-I-21. The average spectrum was extended to 
2.7 GHz and suggested a spectral break with an index of a ~ - 3.0 above 
1.4 GHz. The data collection to determine flux densities of a larger sam- 
ple of pulsars took many years to complete. While numerous observatories 
(Arecibo, Jodrell Bank, Green Bank, Parkes) made observations at frequen- 
cies of 1.4 GHz and below only the Goldstone facility detected pulsars at 
13cm [32]. Observations of the low frequency extension of pulsar spectra 
were carried out in the Soviet Union [18] [19] [94] at frequencies as low as 10 
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MHz. Three pulsars were detected at 8.1 GHz [55]. The suggested existence 
of a spectral break at high radio frequencies [137] was later confirmed [5]. 

A major step forward in the measurement of pulsar flux densities at high 
radio frequencies and hence of pulsar spectra was made by the commission- 
ing of the 100-m Effelsberg radio telescope. Immediately, six pulsars were 
detected at 2.8 cm wavelength [161]. This telescope continued to set records 
of the highest frequencies at which pulsars could be studied by reporting 
detections at 22.7 GHz [10] [11] and finally in the mm-wavelengths [162] [73]. 
A detection of the pulsar PSR 0355+54 at 3 mm wavelength has also been 
achieved with the Pico Veleta telescope [116]. 

The early measurements of flux density of the strongest pulsars had to 
give way to studies of larger samples, if possible with a wide frequency 
flux density coverage. An early compendium of pulsar spectra was given 
for 27 pulsars [141]. Further multi- frequency spectra have been presented 
[7] [142]. In both of these papers the spectral breaks of some pulsars at high 
frequencies were confirmed and in the latter work the frequency evolution of 
the pulse width was noticed (this eventually led to the Radius-to- Frequency 
Mapping concept). Low radio frequency observations [92] confirmed the cut- 
off in pulsar spectra for a number of objects. 

Subsequently, the flux density of a larger sample of pulsars was mea- 
sured at several frequencies [140] [85] [67] and spectra of pulsars were derived 
[93] [153] [106] (Fig. 12). From all these publications the conclusion was that 
the average spectral index is < a > = -18. From our gauss fitted distribu- 
tion (see below) we have found a similar value, < a > = -1.75. It is obvious 
from Figure 26, that the distribution is fairly wide. Some 10% of all pulsars 
require a two power law fit in the high frequency range. A small number of 
pulsars have been recognized with almost fiat spectrum {a > -1.0) [106]. 
In addition pulsar spectra seem to follow the power law down to low fre- 
quencies (a few 10s of MHz) with a few exceptions, where a turn-down is 
observed. 

4-2 Millisecond Pulsars (P < 20 ms) 

Millisecond pulsars were discovered in 1982 [8] as a result of a search in the 
direction of radio sources with very steep spectra. The fiux density of these 
objects is very low. This, combined with the effects of interstellar broaden- 
ing, rendered their detection difficult. Early studies of millisecond pulsars 
[34] [40] suggested that these objects have spectra steeper than slow pulsars. 
Recently high frequency observations of millisecond pulsars were also made 
[66] [107] in spite of their low fiux densities. The spectra of 20 objects were 
studied at lower radio frequencies in a survey of 280 pulsars [85] in the 
northern hemisphere. Several southern millisecond pulsars were also stud- 
ied [153] in a narrow frequency range. A major study of millisecond pulsars 
[74] [167] [75] showed that once a volume-limited sample is considered, many 
of the characteristics of slow and fast pulsars (spectrum, pulse shapes, num- 
ber of sub-pulse components, polarization) are the same. The one distinct 
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difference is found in the luminosity. Slow pulsars arc some 10 times more 
luminous than millisecond pulsars. An investigation of the low frequency 
turn-over of millisecond pulsars [78] [95] revealed that the morphology of 
millisecond pulsars is very similar to that of normal pulsars but with lower 
luminosity. 



5 Polarisation 



Soon after the discovery of pulsars [54] their linear polarization was also 
discovered [89] using the Jodrell Bank Mark I radio telescope. Variations 
of the intensity from pulse to pulse were detected from orthogonal dipoles 
connected to a high speed recorder. The linear polarization was found to 
be surprisingly high even at lower (150 MHz; 408 MHz) radio frequencies. 
Soon it was realized that considerable circular polarization was also present 
in pulsar emission [27] [24]. The early observations with the Parkes telescope 
[125] of the pulsar PSR B0833-45 showed a very high degree of linear po- 
larization of the integrated pulse (in fact nearly 100%) and gave arguments 
for a (magnetized) rotating vector model for pulsars. The observations that 
followed, e.g. [33] [113] , showed that the phase-drift of the linear polarization 
is a common feature in pulsars and, hence, gave support for the rotational 
model. 



5.1 Integrated Pulses 

After these early observations a number of observers embarked on deter- 
mining the detailed polarization characteristics of larger samples of pulsars. 
In all these studies all four Stokes parameters were observed since pulsars 
unlike extra-galactic sources showed high degree of linear and circular po- 
larization. 
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Fig. 13 Early pulsar polarization [96] showing the frequency dependence 
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5.1.1 Normal Pulsars (P > 20 ms) In 1971 the results of observations of 
21 pulsars, at the frequencies of 410 and 1665 MHz were published [96] 
(Fig. 13). Numerous studies of the Crab pulsar in polarization have also 
been made [21] [48] [100] [112]. In 1971 an important result was observed, 
showing that pulsar polarization is constant up to some frequency, after 
which it decreases almost linearly [101]. This was a result not observed in 
any other radio source and required new interpretation. 

Major surveys of polarisation characteristics of larger samples of pul- 
sars started with the access to large radio telescopes like the Effelsberg 
100-m radio telescope [114] [115] [165] [166]. The Arecibo telescope was in- 
volved in polarimetric observations [133] [159]. In the southern skies pulsars 
were studied with the Parkes dish [50] [108] [103] [164] [124] [104]. A multi- 
frequency survey of 300 radio pulsars, at frequencies below 1.6 GHz, has 
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Fig. 14 Polarisation development with frequency of the integrated pulses of nor- 
mal pulsars for the sample of pulsars shown in Fig. 12 [160] 
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Fig. 15 Characteristic polarization behaviour of pulsars. The dark-shaded area 
represents linearly polarized power and the light-shaded area circularly polarized 
power, (a): PSR B0355-h54. (b): PSR 0525-f-21. (c): PSR B1800-21. (d): PSR 
0144-h59. From [155] 



been conducted [47] with the Lovell telescope at JodrcU Bank. The polar- 
ization of a large sample of pulsars was subsequently studied at the highest 
radio frequencies [156] [155] up to the frequency of 32 GHz (Fig. 14). Some 
generalizations about pulsar polarization properties can now be made. 

In Fig. 15 we show the polarization evolution with frequency for four 
characteristic types of pulsars. The pulsar B0355-I-54 (Fig. 15a) begins with 
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low linear polarization percentage at low radio frequencies, then reaches a 
maximum (for one component) in the middle range of frequencies and finally 
falls to low polarization values at the highest frequency so far observed. The 
phase sweep is linear for the highly polarized component but jumps through 
90° between components. Circular polarization is low but evolves in a man- 
ner similar to the linear polarization. Most pulsars evolve in this manner. 
The second evolution sequence is shown for pulsar B0521+21 (Figure 15b). 
Both components are polarized. The degree of polarization falls at high fre- 
quencies and the phase-sweep is S-shaped for both components. In Figure 
15c the polarization-evolution of the pulsar B1800-21 is shown. This has 
the familiar increase and decrease behaviour, but in addition, with a con- 
siderable circular polarization component. A more unusual evolution is seen 
in Figure 15d for the pulsar B0144+59. The linear polarization decreases 
with frequency. However its circular polarization keeps increasing up to the 
highest frequency observed so far. 

5.1.2 Millisecond Pulsars (P < 20 ms) 

The discovery of millisecond pulsars [8] did not lead to immediate studies 
of their polarization. The first reports came in the 90's [150] [117]. This 
was due to the fact that millisecond pulsars are, on average, an order of 
magnitude less luminous than normal pulsars and thus require very sensitive 
polarimeters. A major contribution on millisecond pulsar polarimetry was 
published in 1998 [167] and more recently in 1999 [146]. The polarization 
characteristics of millisecond pulsars are similar to those of normal pulsars, 
namely that the linear polarization falls to high frequencies [167] (Fig. 16). 

5.2 Single Pulses 

Most of the polarization observations presented so far referred to integrated 
pulses. However it is the polarization of the single pulse that tells us about 
constraints that are necessary for the interpretation of their emission mecha- 
nism. Early observations [132] [102] gave strong indications that pulsar radi- 
ation is highly polarized. Soon after a strange behaviour was also detected, 
e.g that orthogonal polarization modes, i.e. abrupt jumps of the position 
angle by 90° in consecutive single pulses are often observed [97] [41]. 

In 1995 a major co-operation project was organized under the auspices 
of the European Pulsar Network. Various telescopes in Europe (Effelsberg, 
Jodrell Bank, Bologna, Westerbork, Torun and Pushchino) were used si- 
multaneously to observe pulsars at a number of frequencies. Recently the 
radio telescopes in Ooty and GMRT (both in India) have joined this net- 
work. Each telescope was optimal at some frequency, so that a very wide 
frequency coverage was achieved. Many of the telescopes have the capability 
to observe the full polarization of pulsars as well. The results of this major 
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Fig. 16 The polarization of millisecond pulsars at I4IO MHz from [167]. Four 
curves are plotted: the total power is shown by the outer curve in arbitrary units. 
The linear polarization is shown by the dashed curve and the circular polarization 
by the dotted curve. The linear polarization position angle is plotted on the lower 
panel in degrees. All curves are plotted against longitude in ms (upper scale) and 
in degrees (lower scale) 



multi-frequency network have shed new light on the wide band performance 
of pulsars as emitters [59] [60] [61] [62] [77]. The fundamental result, that pul- 
sar radio emission is basically broad-band, was confirmed. This is seen in 
Fig. 8 where single pulses observed at four widely spaced frequencies are 
plotted. Observations of full polarization showed even more unusual time 
sequences. While the total intensities correlated rather well, the polarization 
deviations were much higher. In particular the circular polarization (usually 
observed in the conal lobes of a pulse) vary highly from pulse to pulse (Fig. 
17). 



6 Distributions 

More than 1500 pulsars have been discovered to date, facilitating the sta- 
tistical investigation of their distribution in space, period, period derivative 
and in other parameter spaces. These distributions are by now statisti- 
cally stable and reliable, not only because of the large number of stars 
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Fig. 17 Tuio examples of individual pulse pairs of PSR B0329+54, observed in 
full polarization. The dark-shaded area represents linearly polarized power and the 
light-shaded area circularly polarized power. From [59] 

involved in the sample but because they cover a large portion (~1%) of 
the pulsar population which is believed to be of the order of 10^ pul- 
sars in our Galaxy. The figures presented below were produced using the 
data for 1300 pulsars available in the recently released ATNF pulsar cat- 
alogue (www.atnf.csiro.au/research/pulsar/psrcat/) [105]. A few pulsars in 
the ATNF pulsar catalogue are X-ray pulsars. Among them there are some 
with long period. They have been included as their properties arc very sim- 
ilar to normal radio pulsars. 

6.1 Spatial distribution 

The distribution of pulsars in galactic coordinates is shown in Fig. 18 in 
Hammer- Aitoff projection. It is obvious that pulsars are strongly grouped 
along the galactic plane. Millisecond pulsars (many of which are also binary) 
are more isotropically distributed. This effect is due to their inherent weaker 
emission, which allows the detection of, primarily, nearby objects. 

6.2 Dispersion Measure vs. Galactic Latitude distribution 

The Dispersion Measure (DM) of pulsars depicts the electron content on the 
line of sight, between the star and the observer, J ne{l) x dl, where ne{l) 
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Fig. 18 The distribution of 1395 pulsars in Galactic Coordinates. Normal pulsars 
are depicted by dots. Binary pulsars have a large circle around them. Pulsars with p 
< 20 msec are depicted by stars. Binary millisecond pulsars are shown by encircled 
stars. Figure courtesy of B. Klein 

is the electron density and 1 is the distance. If the average electron density 
is known, then the DM is a direct measure of the pulsar distance. From 
Fig. 19 it is obvious (a) that pulsars are tightly clustered along the galactic 
plane (galactic latitude, b = 0°) and (b) that the highest DMs are found on 
this plane. This is due to the electron density distribution in the Galaxy, 
which has been modelled by several researchers [91] [148]. According to these 
models, the electron density peaks at the Galactic centre and falls off with 
a scale height of about 70 parsecs away from the plane. 

A few pulsars deviate from the smooth 1/e distribution. These objects 
are known to be located behind dense HII regions, c.g the Gum Nebula. 



6.3 Period distribution 

Pulsar periods range from 1.5 ms to about 8.5 s (Fig. 20). This range of 
values for a physical parameter characterizing one species of objects is too 
wide to be explained by a common origin. It is widely accepted that ms 
pulsars are recycled pulsars, spun up by accretion processes, during which 
they accumulate mass and obtain extra spin (from this mass) from a binary 
component. This is confirmed by the data of Fig. 18, which show that most 
millisecond pulsars are in binary systems. 

The histogram of pulsar periods shows a distinct bimodal distribution. 
The median of the period distribution of normal pulsars is about 0.65 s, 
whereas for millisecond pulsars it is 0.0043 s (4.3 ms). There is a character- 
istic lack of pulsars with period around 20 ms. 
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Fig. 19 The distribution of the Dispersion Measure of 1282 pulsars as a function 
of their Galactic Latitude 



6.4 Period derivative (P) distribution 

The P distribution of pulsars (not shown, as it is directly connected to Figs. 
20 and 21) varies between 10"'^^ ss^^ to 10^^° ss^^. Millisecond pulsars 
exhibit a much slower decay (lengthening) of their period. One of the fastest 
decaying period pulsars is the Crab pulsar, whose period slows down by 36 
ns per day. In general the P distribution of pulsars is very similar to the 
period distribution. 

There is a noticeable bimodal distribution. The median of the P distri- 
bution of normal pulsars is 2.6 x 10^^'^ ss~^ , whereas for millisecond pulsars 
it is of the order of 2 x 10^^° ss^^ . Five orders of magnitude slower period 
decay. 



6.5 Period - Period derivative (P) distribution 

One of the most important graphical distributions of pulsars is their period 
(P) - period derivative (P) distribution (Fig. 21). As expected from Figures 
20 and 22 it shows very clearly the characteristic clustering of normal pul- 
sars (large P, large P) and of millisecond pulsars (small P, small P). The 
vast majority of normal pulsars are isolated single stars. On the contrary, 
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Fig. 20 The period distribution of 1300 pulsars 



the majority of millisecond pulsars are members of binary star systems. Bi- 
nary pulsars located between the two clusters will slowly drift toward the 
millisecond cluster in less than 10® years. 

Assuming that the decay of pulsar periods is due to their dipole radia- 
tion, their characteristic age can be calculated from the very simple expres- 
sion Tchar — \^ years. The "dash-dot" lines in Fig. 21 correspond to lines 
of constant age. The characteristic age of normal pulsars is of the order of 
10* years, whereas the age of millisecond pulsars is slightly above 10^ years. 
The Crab pulsar is the isolated pulsar closest to the 10'^ years line. 

Following classical electrodynamics theory, the surface magnetic field 
of pulsars is given by the expression Bq = 3.3 x 10^^ x \/ Pp gauss. The 
"dash" lines in Fig. 21 correspond to lines of constant magnetic field. It 
is immediately noted that normal pulsars have a surface magnetic field of 
about 10^^ gauss, whereas the surface magnetic field of millisecond pulsars 
is much lower, of the order of 10^ gauss. 

Finally it should be mentioned that the absence of pulsars in the lower 
right corner of the diagram is due to the existence of a "death line" , owing 
to the gradual decaying of the induced electrical potential of pulsars. Slow 
pulsars with low magnetic field cannot develop a large enough potential 
above their magnetic poles for discharges (and therefore radiation) to take 
place. The absence of millisecond pulsars below about 10^" years indicates 
that their age cannot exceed the Hubble time (age of the Universe). 
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Fig. 21 The P - P distribution of 1395 pulsars. Binary pulsars are depicted by 
encircled dots 

6.6 Age distribution 

The characteristic age {Tchar ~ \ ^) distribution of pulsars (Fig. 22) shows 
the expected bimodal distribution attributed to the different ages of normal 
pulsars and millisecond pulsars. The median age of normal pulsars is 4.7 x 
10^ years, whereas the age of millisecond pulsars is of the order of 5 x 10^ 
years. 



6. 7 Surface Magnetic Field distribution 

The surface mag nctic field (Bq = 3.3 x 10^^ x ^PP) distribution of normal 
pulsars is tightly peaked at 1.3 x 10^^ gauss, whereas for millisecond pulsars 
is much lower, of the order of 4 x 10* gauss(Fig. 23). In a few stars the surface 
magnetic field is larger than 10^"' gauss. This is the surface magnetic field 
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Fig. 22 The Age distribution of 1194 pulsars 



of a new class of objects, called magnetars, first detected through their X- 
ray emission [69]. Magnetars are known to be neutron stars, which do not 
always emit at radio wavelengths. 



6.8 Luminosity distribution 

The luminosity distribution at 400 MHz of 612 pulsars is depicted in Fig. 24. 
The mean of the gauss fitted distribution is 115 mJy kpc^. The luminosity 
is calculated from the 400 MHz flux density, assuming that the Dispersion 
Measure is a true measure of the distance of each pulsar. This may lead to 
an over estimate of the luminosity of pulsars located behind HII regions. 



6.9 Spectral Index distribution 

The spectral index distribution between 400 MHz and 1400 MHz of 285 
pulsars is shown in Fig. 25. It is a rather wide distribution. The bulk of 
pulsars demonstrate a spectral index between -3 and 0. The mean of the 
gauss fitted distribution is -1.75 ± 0.1. Pulsars with flat spectral indices are 
the ones which should be investigated at high frequencies. 
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Fig. 23 The surface magnetic field distribution of 1194 pulsars 
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and 1400 MHz were used 
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